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Abstract 

We study chaotic eigenfunctions in wedge-shaped and rectangular regions using a generalization 
of Berry's conjecture. An expression for the two-point correlation function is derived and verified 
numerically. 
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I. INTRODUCTION 



According to Berry's conjecturej^ chaotic eigenfunctions behave locally like random su- 



perpositions of plane waves with wavevector k, where k = (l/ti)\/E — V. This description 
is consistent with random matrix theory. The Gaussian random wave model does not ac- 
count for the localization properties of eigenfunctions, such as scarringi and weak quantum 
ergodicity.i In particular, for quantum billiards with Dirichlet boundary conditions, the 
wavefunction must vanish at the boundary, and thus no longer looks random in its vicinity. 
For straight boundaries, the problem can be solved by reducing the set of available plane 
waves to those which are antisymmetric with respect to reflection across the wall, thereby 
ensuring a zero value along that boundary; see Berry.i If, for instance the boundary is given 
by the line y = 0, the chaotic wave function will be composed of a seemingly random combi- 
nation of plane waves of the form sin(k y y)cos(k x x + 4>) where k^ + k 2 , = k 2 and is a random 
phase shift. Bies and Helleri discuss similar boundary effects in soft potentials. 

Some criteria for unconfmed random waves appear to be well satisfied by eigenfunctions 
of chaotic billiards. McDonald and Kaufmani and McDonaldl checked for random nodal 
patterns, Gaussian statistics and Bessel function correlations. In their original work on 
the stadium, these authors emphasized the qualitative appearance of nodal lines in the 
eigenfunctions. The nodal lines appear to wander randomly throughout the billiard, which 
is indicative of an isotropic distribution of local wave vectors. Furthermore, they found the 
wavefunction statistics to be Gaussian. However, the correlation function is not a Bessel 
function for a general billiard. 

We first generalize Berry's boundary-adapted form of the Gaussian random wave model 
to wedge-shaped regions. The case of a 90° angle is immediately solved: one just antisym- 
metrizes with respect to reflections in both x and y. The simplest non-trivial case consists 
of a wedge with an opening angle of 60°. We focussed on this example in what follows, 
although our method extends to all opening angles of 7i/n radians, where n is a positive 
integer. We first determine the right boundary-adapted plane wave basis. The two-point 
correlation function can be used to compare our model to numerically generated chaotic 
eigenstates. For Gaussian random waves in free space, the later is known to give 

MxMx')> = Jo(|x-x'|) (1) 

where the average is taken over a set of eigenstates. As we shall see, the two point correlation 



function for 60° wedge again reduces to a sum of Bessel functions. Our numerically generated 
ensemble of wavefunctions is compared to this form of the correlation function, including 
the interference effects of two or more Bessel functions. 

Our second type of billiard consists of two semi-infinite parallel lines at x > 0, y = ±a/2, 
which are connected by a perpendicular wall at x — 0. The wavefunction can be made to 
vanish on the back wall by antisymmetrizing it with respect to the y-axis. The parallel 
walls are not so readily handled. The antisymmetrization procedure extends over an infinite 
periòdic array of such lines a distance a apart. The two-point correlation function will 
be composed of an infinite sum of Bessel functions. Since the later diminsh rapidly with 
distance, only the nearest neighbours contribute, allowing us to compare our formula with 
its numerically determined value. 

II. THE TWO-POINT CORRELATION FUNCTION 
A. The wedge 

We first look for the properly antisymmetrized wavefunctions. Let 

V>(x) = / d6a e e- lke * +l5e (2) 

denote a random sum of plane waves. Here kg = A;cos#x + ksin8y, the ag are independent 
Gaussian distributed random variables and the 5o are independent uniformly distributed 
random phase shifts. Let denote the desired boundary-adapted version of ipfe). Since 
^(x) vanishes on the boundary of the wedge, it can be extended to a fictitious wavefunction 
living in the whole plane by reflecting antisymmetrically across either boundary. To obtain 
ip(x), we let Ri denote reflections with respect to one of the lines bounding the wedge and 
R 2 reflections with respect to the other line. We would like to antisymmetrize ^(x) with 
respect to both of these reflections. The projection operator (1/4) (1 — Ri)(l — R2) does not 
work, however, because R\ and R 2 do not commute. Instead we must resort to the following 
group-theoretical construction. Let R 3 denote the reflection with respect to the line that 
passes through the vèrtex of the wedge and meets the two edges of wedge at a 60° angle. 
The product R\R 2 of the two reflections R\ and R 2 is a rotation through 120° centered on 
the vèrtex, and R 2 R\ is the opposite rotation. Thus, the set 1, R±, R 2 , R3, R\R 2) R 2 R\ forms 
a representat ion of the dihedral group C3. It is now easy to see that V»(x) transforms under 
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the one-dimensional representation of C3 that assigns a character of -1 to the reflections 
and 1 to the rotations. Put another way, the problem of going from ip(x) to ip{x) is that 
of projecting onto this irreducible representation of C3. Formally, this transformation is 
expressed as! 

^( x ) = "i E X(A)A^). (3) 
V A6 c 3 

Here, A denotes the group elements 1, Ri, R 2 , R3, R1R2 and R2R1, and the character is given 
by x(Ri) = x{Rj2) — XÍ.R3) — ~ 1 an d = XÍRiRq) = xi^Ri) = 1- The normalization 
factor of 1/ is chosen to ensure that the two-point correlation function gives the free-space 
result far away from the wall. If we then restrict ^(x) to the wedge-shaped region, we obtain 
the required boundary-adapted sum of plane waves. 

Now that we have an explicit expression for ip(x) we can use it to calculate the two-point 
correlation function C(x, x + r) = ( , 0(x)'0(x + r)). The correlation function is a sum of 
thirty-six terms of the form 

( X (A)A^)x(B)B^ + v)) = X (A)x(B) J d 6e~^ ^- B ^ Br (4) 

= x(A)x{B) J d6é kpcos{e - do) (5) 
= X (A)x(B)J (kp) (6) 

Here p = |v4x — B(pc + r)\, that is, the distance from x + r to B~ l Ax. The correlation 
function becomes 

C(x,x + r) = ^ E x(A)x(B)J (kp B - ÍA ) (7) 

A,BeC 3 

= E x(C)J (kp c -i) (8) 

C&Cz 

= J (kr) - J (kp Rl ) - J (kp R2 ) - Jo(kp R3 ) + J (kp Rl R 2 ) + J (kp R2Rl ) (9) 

In the second line, we replaced B by AC . The sum over A then becomes trivial, since 
x(A)x(B) = x(A)x(A)x{C) = x(C). In the last line pa denotes the distance from x + r to 
Ax. 

This formula predicts that the two-point correlation function should display the inter- 
ference of Bessel functions centered at all the points given by applying the group transfor- 
mations on x. Far from the boundary, Eq |f] reduces to the single Bessel function Jo(kr), 
the value that it would have if no boundary were present. However, the interference is 
pronounced near the boundary, and more especially near the vèrtex of the wedge. 



B. the rectangular corridor 



In the case of the rectangular corridor, the back wall and the semi-infinite walls need 
to be handled separately. The wavefunction at x in the x < half-plane is obtained by 
reflecting about the y-axis. To extend the wavefunction in the y direction, we reflect it with 
respect to the wall at y — a/2. By repeating this process for all y — (1 + 2m)a/2 where 
m = {..., —1, 0, 1, ...}, the whole plane is tiled with positive or negative copies of 

Mathematically, the process by which ^(x) is obtained is as follows. The initial billiard is 
placed on an infinitely long cylinder of circumference 2an. Our initial x-axis is now parallel 
to the axis of the cylinder, whereas the y-axis is wrapped around the cylinder and truncated 
at y — (1 ± 2n)a/2. The wavefunction is first reflected across the y-axis. Therefore the 
antisymmetrization in the x-direction transforms as a one-dimensional representation of the 
reflection group Z 2 . Rotations by múltiples of 2a should give back the wavefunction, while 
reflections about any of the lines parallel to the x-axis give its negative version. We thus 
obtain a representation of the dihedral group C n . It is now clear that on the cylinder, ip(x) 
transforms under the one-dimensional representation of C n <g> Z 2 that assigns a character of 
— 1 for reflections and 1 for rotations. This time the transformation is given by as 

1 

" Aec„®z 2 

As in the case of the wedge, the normalization factor is found by requiring that the correlation 
function give Jo(kp) far away from the wall. 
The correlation function is given by: 

C(x,x + r) = lim^oo- 1 - Yl x(A) X (B)J (kp B -i A ) (11) 

4n A,B£C n ®Z 2 

= lim^oo Yl x(C)Jo(k Pc -i). (12) 
cec n ®z 2 



^(x) = -= £ X(A)A^). (10) 

V4n 4p/7„68Zo 



However, J (kp c -i) — > \j \/kp c -\. Hence, only the Bessel functions centered at nearby 
points contribute, and the correlation function remains finite as we take the n — > oo limit. 



III. NUMERICAL RESULTS 



Here we apply the theory to check further the properties of random waves in billiard 
systems. There are some obvious and some more subtle modifications of random wave 
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behavior known in closed chaotic billiards. The wavefunction must vanish on the boundary, 
and scarring affects some states in a nonrandom way.i The theory given here suggests new 
correlations which ought to be checked in chaotic systems. As we incorporate more about 
a specific billiard geometry into our correlation functions, we are probing the properties of 
waves which are "as random as possible" within the constraints, such as a wedge boundary 
being present. In checking the numèrics of a billiard such as the cone or stadium, we are 
really asking whether the eigenstates are truly random but for the constraint of interest, say 
a wedge boundary. Naturally this cannot be strictly true, since there are more constraints 
that we haven't included. At the end of this strategy comes a tautology: if, in a closed 
billiard we incorporate all the constraints, we have only the eigenstates left, as the only 
waves which are consistent with all the constraints! 

A. The wedge 

To validate our expressions for ip(x) and C(x, x + r), we generated an ensemble of 500 
eigenstates near k = 200 for a 60° wedge which has been closed off by a semicircle; see 
Fig 1. The eigenstates were found using the Boundary Integral Method. Our cone-shaped 
billiard has a circle diameter equal to 1, which is about 32 wavelengths across. The Poincare 
section indicates that this billiard is chaotic. In Fig. 2, we show C(x, x + r) for (a) x on the 
symmetry line through the center of the billiard, and (b) for x placed within a wavelength 
of one of the edges. Both Fig. 2(a) and 2(b) display the expected interference, but the later 
is especially pronounced in Fig. 2(b). Note that since x is on the symmetry line, only even 
states contribute to Fig. 2(a). 

For comparison, Fig. 3 displays the theoretical predictions for the two-point correlation 
functions of Fig. 2. Both plots are in good agreement with their numerical equivalent. In 
particular, the interference pattern of Fig. 2b is reproduced in Fig. 3b. However, the angular 
nodes are much more pronounced in the theoretical plots, suggesting that an unaccounted- 
for smoothing process is at work in the numerical experiment. We compared the error 
quantitat i vely using 

/Gpr|C num (x,x + r) -C th (x,x + r)| 2 

Jrf 2 rC t h(x,x + r)2 ' [ > 

We obtained 0.52 for Fig. 2a and 0.15 for Fig. 2b. The discrepancy in the predicted and 

experimental vàlues is expected due to the difference in the angular nodes, and appears to 
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be more pronounced far away from the billiard walls. 

In Fig. 4, we take the angular average of the correlation functions of Fig. 2. As was found 
by Li and Robnikjïl Fig. 4a resembles the Bessel function J (kr). This should be the case 
whenever boundary effects contribute equally from all sides. Fig. 4b also approximates a 
Bessel function because, as can be seen from Fig. 2b, the boundary effects affect only a small 
range of angles and furthermore are independent of radius. 

B. The rectangular corridor 

Once again, we took an ensemble of 500 eigenstates near k = 200, this time for a quarter 
stadium with a circular radius of R = 0.6 and straight length l = 1.2 (about 20 wavelengths 
wide); see Fig. 5. The numerical correlation function is shown in Fig. 6 for x (a) on the 
symmetry line through the center of the billiard, (b) within 1 À of the top wall, (c) within 
1 A of a còrner and (d) within 1 À of the back wall. As expected, only the Bessel functions 
which are less than a few wavelengths away contribute to C(x, x + r). 

The numerical correlation function was compared to the theoretical prediction using 
Eq. [H| Here only the first few terms were kept in C t h(x, x + r). We obtain 0.55 for Fig. 6a, 
0.23 for Fig. 6b, 0.03 for Fig. 6c and 0.24 for Fig. 6d, in agreement with our results for the 
cone billiard. 

IV. CONCLUSION 

We have successfully extended the boundary-adapted Gaussian wave model to a wedge- 
shaped region with opening angle 60°. Our technique is readily generalized to any opening 
angle of the form ir/n for n integer; one merely replaces the group elements of C3 with those 
of the dihedral group C n . We also solved the case of the semi-infmite rectangular corridor. 
It is our hope that these result will stimulate more work on boundary effects in arbitrarily 
shaped billiards. 
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FIG. 1. The cone billiard. 




FIG. 2(a). Experimental correlation function for x = (0.3, 0), on the symmetry axis through the center of 
the cone. The grid is a square of side 0.129, or about 8 wavelengths, centered on x. 




FIG. 2(b). Same as Fig. 2(a) for x = (0.3,0.153), near the upper straight edge of the cone. 
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FIG. 3(a). Theoretical correlation function for x = (0.3,0) (same grid as in Fig. 2(a)). 




FIG. 3(b). Same as Fig. 3(a) for x = (0.3,0.153) (same grid as in Fig. 2(b)). 




0.14 



FIG. 4(a). Angular average of C(x + r,x) vs r for x = (0.3,0) (same grid as in Fig. 2(a)). Sòlid line, 

Bessel function J (fcr); x, numcrical data. 
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FIG. 4(b). Same as Fig. 4(a) for x = (0.3,0.153) (same grid as in Fig. 2(b)). 
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FIG. 5. The stadium billiard. 




FIG. 6(a). Experimental correlation function for x = (0.3,0.3), the mid-point of the square part of the 
stadium. The grid is a square of side 0.140, or about 8 wavelengths, centered on x. 
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FIG. 6(b). Same as Fig. 6(a) for x = (0.0223,0.2), near the back wall. 




FIG. 7(a). Theoretical correlation function for x = (0.3,0.3) (same grid as in Fig. 6(a)). 




FIG. 7(d). Same as Fig. 7(a) for x = (0.2,0.0255) (same grid as Fig. 6(d)). 
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